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A Lorentz covariant kinetic equation for bound states and their constituents is presented 
and solved exactly in closed form. It describes in a unified way dynamical formation and 
£SJ \ dissociation of states such as quarkonia and (anti)-deuterons in the excited medium formed 

; I ■ with a high-energy heavy-ion collision. 
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The possibility to describe hadronic bound states in strongly interacting matter provides an impor- 
tant tool for understanding the dynamics of the medium produced in high-energy heavy-ion collisions. 
The main reason lies in the fact that bound states can be formed and destroyed during specific stages 
of evolution of the medium, with a delicate interplay among various elements: the binding energy and 
structure of the object under consideration, the medium nature and characteristics which change as 
function of time and the abundances and phase-space correlations of the bound state constituents. 

There are essentially two important kinds of bound states which deserve attention: heavy quarkonia 
(J/?/>, T and excited states) and (anti-)deuterons (d, d). Their properties are quite different from one 
\ another and for this reason they probe different stages of a nuclear collision. The former are tightly 
bound, small objects, in general characterized by a hard scale and subject to confinement. They probe 
the early stages of a collision, in particular a possible quark-gluon plasma (QGP). The latter are very 
loosely bound, large objects which can be formed and survive dissociation only at the very latest 
q ■ moments of a collision. They probe the freeze-out stage. 

Despite the aforementioned differences, the essential physics which regulates the time evolution in 
an interacting medium is very similar. For both bound states it is possible that they are formed in 
the evolving medium and are subsequently dissociated. How this happens is regulated by transition 
^ ■ probabilities and by the medium constituents' phase-space density. It is therefore desirable to have a 
unified description of these phenomena. 

In the following we discuss a general kinetic equation, appropriate for both cases, and we present in 
closed form the exact solution, examining in detail specific features of the two cases under study. The 
solution is constructed for two different ways of specifying the initial condition. Finally, we discuss a 
simple analytical example, leaving the more realistic cases for future numerical calculations. 

We begin by specifying the processes B(p) + n(k) «-> c\(qi) + 02(92) of dissociation and formation 
of a bound state B into and from a pair of constituents C\C2 by means of a fourth degree of freedom 
labelled n. The corresponding momenta are indicated in parentheses. As mentioned earlier, the bound 
state B can be a quarkonium or a(n) (anti-)deuteron, while the respective constituents are qq pairs of 
heavy quarks or proton-neutron pn pairs. In particular, the dominant reactions to be considered are 
cc <-> J/ipg(ir), bb <-> Tg(n), pn <-> ^(7), pfi <-> dirfa). The amplitude for the process is illustrated in 
Fig.ffl 

Borrowing concepts from the treatment of the non-linear Boltzmann equation ,1] we postulate that 
the kinetic equation describing the time evolution of the phase-space density Tb of the bound state 
is then 

P^d^ Fb(p, x) = C F (p, x) - C D (p, x) Fb{v, x) , (1) 
consisting of a drift term (the differential operator p^dn) and a collision term, the latter made of a 
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FIG. 1: Diagram illustrating the amplitude for the formation of a bound state. The dissociation amplitude is 
obtained with detailed balance. Possible reactions are cc <-> J/tpg(7r), bb <-» Tg(7r), «-> ^(7), <-> dir(j) 



formation part 

C F (p,x) =J2 d $(p,k,q 1 ,q 2 ) W C1 (s,t) (2) 

n 

X fc 1 ( ( ll^ x )fc2(Q2,x) 

and a dissociation part 

C D (p,x) = y / d$(p,k,q 1 ,q 2 ) WBn-,c lC2 f n (k,x) . (3) 

In the above / Cl , / C2 and / n are the phase-space distributions of the degrees of freedom participating in 
the collisions, depending implicitly on collision energy and centrality, while the phase-space integration 
measure is 

.,. , . 1 d 3 k d 3 qi d 3 q 2 , ., 

d*( P ,k, qi , q2 ) = - {27r)32Ek {27t)32Ei {27t)32E2 (4) 

x (2T,) A 5 A (p+k- qi -q 2 ). 

As evident from the integration measure, all particles are assumed to be on their mass shell. The 
transition probability satisfies detailed balance as 

W ClC2 ^ B n{s,t) = WBn->cic 2 (s,t) (5) 

and a sum over all particle species n constituting the medium and contributing to dissociation and 
formation is indicated. Clearly, which degrees of freedom are active depends on the evolution time. 
The arguments of the transition probability are the Mandelstam variables s = (p + k) 2 = (qi + q 2 ) 2 
and t = (gi — p) 2 = (k — q 2 ) 2 . In the equations above we have neglected Bose-Einstein or Fermi-Dirac 
statistics, i.e. Bose enhancement and Pauli blocking. 

All the dynamics is well defined once one specifies the distributions of the medium constituents 
and those of the bound state constituents, together with the transition probabilities W ClC2 ^Bn- The 
kinetic equation can then be solved in a closed form, which is convenient to present in two different 
ways, depending on how the initial condition is given, i.e. either at some initial global time to or 
at some initial proper time To- The form of eq. Q is very general and in the following we make no 
reference to the detailed structure of the collision term. In other words the functions Cf and Cb are 
assumed to be known. 

When the initial condition is specified at some initial global time to to be J- B (p,r), then the kinetic 
equation is conveniently re-written as 

d d\ 

-Q t + v ■ — j F B (p, f, t) = \ F (p, r, t) (6) 

-\ D (p,r,t) F B (p,r,t) , 
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where v = p/E p and Xf,d = Cfd/E p . The solution to this last equation can be obtained with some 
simple steps. Neglecting the collision term, it is trivial to see that the solution is the scaled function 



T* B {p,r,t)=J^ B (p,p(t-to)) , 



(7) 



where p{t) = r — vt. This is the case of free streaming. The case with dissociation but without 
formation was studied in |2J to address the problem of charmonium suppression. There the solution 
was given as 



(8) 



x exp 



to 



dt'X D (p,p(t-t'),t') 



Otherwise, with the only formation term and without dissociation, the solution is obtained by direct 
integration as 

(9) 



F$ im (P, r, t) = F* B (p, f, t)+f dt' X F {p, p(t-t'),t') 



It is interesting to notice the structure of the scaled argument p{t — t') = r — v (t — t') , ensuring that 
any function of it is annihilated by the drift operator. This observation, although trivial, is crucial for 
the complete solution. We now combine the previous partial solutions and seek for the general case a 
solution of the form 



T* B (p,r,t)+[ dt'a(p,p(t-t'),t' 

f dt' X D (p,p(t-t'),t') 



x exp 



(10) 



where a is a function to be determined. With this choice, substituting in eq. © one finds right away 
that the needed function is 



a(p, f, t) = X F (p, f, t) exp 



t 



dt'X D (p,p{t-t'),t) 



(11) 



Then, the general solution of the kinetic equation as function of global time can be obtained by 
substituting this last result into the trial form of eq. (fTT|) . The final result is then 



T B {p,r,t) = I T B (p,f,t) exp 



to 



dt' X D (p,p(t-t'),t') 



(12) 



+ dt' X F (p, p{t - t'),t') exp - jf dt" X D (p, pit - t"),t") | . 



'to 

By immediate inspection one finds that the initial condition is satisfied at t = to and that the two 
previous partial solutions J- B ss in eq. (jHJ) and Jg™ in eq. Q are obtained in the respective limiting 
cases X F = and Xd = 0. We will comment on the obtained result later on, after the second derivation 
of the solution 

Before we continue note that integrating over the position variables gives the invariant spectrum 



E, 



dN B 



d 3 i 



p d 3 p J (2vr) 3 
of the bound state at a given final global time tf. 



F B (p,r,t f ) 



(13) 



4 



If the initial condition is specified at some initial proper time To it is convenient to use momentum 
rapidity y = 1/2 log ((-E + p z )/(E — p z )} and space-time rapidity i] = 1/2 log ((i + z)/(t — z)) as 
longitudinal variables. The initial condition is then J-^(y,p±, r), f±). It is then useful to re-write the 
drift operator as p M <9 M = V = Vn + T>± where 



Pii = m± 



cosh(y - r,) A + I sin h( y - 77) A 



(14) 



V ^-WL (15) 

As in the previous case, let us first consider the situation of free streaming, without a collision term. 
Neglecting the transverse coordinate dependence for a moment, the solution of the free streaming 
equation with the only operator Py is any function of the variable 

f = — sinh(y - rj) , (16) 
to 

It is therefore convenient to define the quantity 

Y(0=y-log^ + VTTe) , (17) 
which reduces to rj when t = tq. In this way one obtains the longitudinal free streaming solution 

F* B \\(y,v,T)=?t ll {y,Y(0) . (18) 

Concerning the transverse part we look for a scaling solution of the type 

H±(r±, r) = F° B± (f ± - v ± f) , (19) 
with v± = p±/m± and f a function to be determined. The kinetic equation becomes 

vn± = g|- f B i ■ <-"±) vt t + v±n± (20) 

which is equivalent, in the non-trivial case V>^T^^ ^ 0, to the equation f = m±. It is simple to 
see that the solution is 

TK - T) -f n "Wrkp? (21) 

with the correct boundary condition f = when t = tq. The full free streaming solution is obtained 
by putting together the separate results for the longitudinal and transverse parts, yielding 

n(v,P±,V,f±,r) = ^ B {y,p ± ,Y(0,p±(^,r)) (22) 

being p±(£,r) = f± — v±t(£,t). We now consider the more complicated case with dissociation, still 
neglecting formation. Analogously to the previous case, we look for a solution of the form 

^B SS (y,P±,V,r±,r) = T* B {y,P±,V,r±,r) (23) 

x exp [-E(y,p ± ,f],r ± ,T)] , 

where the exponent E is determined by the equation 

VE(y,p±,r],r±,T) = C D (y,p±,r],r±,T) . (24) 
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The latter can be solved by direct integration, taking care of the fact that the integration measure 
should be the same as that of eq. (|21l) . ensuring the correct scaling with the variable £. The solution 
is then 

E{y,p±,v,r±,T) = dr = ^ (25) 

J T0 m ± v 1 + (£Vt') 2 

With this last result, especially in its structure concerning the integration measure, it is now possible 
to obtain the full solution. Looking at the earlier solution given with eq. (|12|) in the previous case of 
initialization at fixed global time, one can repeat the arguments leading to it and obtain, after defining 
the rates Ap,D = Cf.dA^x for dimensional clarity ([Ajr£>] = fm _1 ), the general solution of the kinetic 
equation as 



FB(y,P±,ri,r±,T) = |j^(y,pj_,r/,fj_,r) exp 

f T , , ^F(y,P±,Y^T /r'),p ± (C,T')) 
+ I dr exp 



j A D {y,p±,Y^T /r'),p ± (tr')) 

Vi + (£Vr') 2 

_„ A D {y,P^Y(Cr /r"),p ± (tr")) 
y/1 + (£t /t") 2 



(26) 



In this case the initial condition is satisfied at r = To and the partial solutions J-g ss and J r ^ vm are 
again recovered in the respective limiting cases A^ = and Ad = 0. 

The structure of the solution is rich of details on the dynamics of dissociation and formation. The 
first part of the r.h.s. of eq. ((2*6^1 . but see also eq. (|12[). gives the final phase-space density of those 
bound states which were initially present at the start of the evolution at r = tq. This term is especially 
important for the study of quarkonia. These are initially produced by hard collisions when the two 
nuclei first interact and are subsequently dissociated by the formed medium. On the other hand, at 
the high collision energies considered, (anti-)deuterons are absent at any early stage of a heavy-ion 
collision, therefore this term can be neglected. The second piece of the solution describes formation 
of bound states from tq up to t' and their subsequent dissociation from r' until r. For quarkonia 
it describes formation in a QGP, process that might be significant when the phase-space occupation 
of heavy quarks is large enough. This might happen at high collision energies, as recently suggested 
0, 0]. For (anti-)deuterons this second term is the relevant one, describing coalescence of pn pairs. 
This process is relevant only at low enough medium densities (small Ao), i.e. at the freeze-out stage. 
A more sophisticated treatment of this problem for deuterons at low and intermediate collision energies 
was developed in Here, since we consider the high-energy collision regime, we can safely neglect 
many-body correlations. Hence the problem is significantly simpler. 

An interesting application of eq. (|26|) for the study of quarkonia is the comprehensive analysis 
of rapidity [j| and transverse momentum [(| dependencies of medium effects on the final spectrum. 
Regarding (anti-)deuterons, collective flow patterns emerge naturally as they are built-in both in 
formation and dissociation rates, providing a useful probe of the freeze-out properties of the produced 
medium jjfl. Details of how different effects play a role both for quarkonia and for (anti-)deuterons 
can only be addressed with sufficiently realistic numerical studies. In fact, a detailed knowledge of the 
medium constituents' and of the transition probabilities is required. In the following we will limit our 
discussion to a simple example to clarify the physical content of the result contained in eq. (|26j) . 

Before doing so we recall that the invariant spectrum of bound states can be computed by means 
of the Cooper- Frye formula ja|. Given a hyper-surface ry = Tf(rj,f±), one has 



dN B 
dy d 2 pj 



(2vr) 3 



Fb (y,P±,V,r±,T f (rf,r±)) 



(27) 
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where Tb is given by eq. (|26|) and 

dc^p 11 = d 2 f±dr] Tf m± cosh(y — 77) (28) 

- dr f dr f ■ w s 

-T f p±--^-- m ± —smh(y - rf) 

is the invariant volume integration measure. 

To discuss the physical content of the solution given with eq. (|26|) . or equivalently with eq. (|12[) . we 
focus on a single aspect of it. We refrain from discussing the part describing dissociation of initially 
present bound states. This was extensively done in 0| for charmonium and only specific models of the 
medium can give qualitatively new answers. On the other hand the part describing formation deserves 
more attention. As mentioned earlier it can potentially account for quarkonium formation in a QGP, 
as suggested in 0] and become the dominant mechanism of production at high enough energies (RHIC, 
LHC). Moreover, the second term should describe formation of (anti-)deuterons in the late stages of 
the collision. While much work has been done in this domain, efforts were limited to applications of 
the coalescence model at freeze-out (See 0,0] for some recent results). Only recently 0] a new 
attempt was made to address this problem. 

Therefore, focussing attention on bound state formation in the medium, we neglect the first term 
in eq. (|26|) either because dissociation destroys all the initially present bound states or because there 
are no bound states to start with. With significant simplifications of the solution in order to follow 
the essential physics, we neglect all dependencies on y,p±,r] and f±. The remaining solution can be 
written as 



N b {t)= / dr' K f {t') exp - / dr" A d (t 

Jtq I Jt' 



(29) 



This situation is equivalent to having solved the first order ordinary differential equation dNs/dr = 
Af — Ad Nb with initial condition Ng = 0. We then assume that the time dependence of the formation 
and dissociation rates are 

A f (t) = PfNA and A D (r) = ^, (30) 
r r 

where Pfd are formation and dissociation probabilities, N Cl 2 indicates the number of particles which 
can form the bound state (either pn or qq pairs) and Njy is the number of particles that can dissociate 
the bound state. With the latter assumption one can readily perform the time integrations in eq. ()29|) 
obtaining the final number of bound states at tj as 



N 



f _ P F N C1 N, 



(■■2 



B 



P D N D 



l-(r /r f ) P - N - (31) 



when PdNd 1. This last results indicates that the final yield is proportional to the number of 
particles which can form the bound state and inversely proportional to the number of particles that 
can dissociate it. With the help of numerical simulations, the consequences of this scaling law for 
deuterons were discussed previously in [13j. Moreover, since the number of dissociating particles is 
proportional to the volume of the medium, either at hadronization for quarkonia or at freeze-out 
for (anti-)deuterons, we recover the well known coalescence formula here newly obtained with 
dynamical considerations. The above simplified derivation gives some insight in the coupled dynamics 
of dissociation and formation. In particular, as discussed in [lj, LLD] , the experimental observation 
of an inverse proportionality of the number of quarkonia with the total number of produced hadrons, 
might give support to this new mechanism of production in a QGP. 

We can now conclude by summarizing the results presented above. The aim of this paper has 
been twofold: to present a unified theoretical description of bound state formation in the strongly 
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interacting medium, formed with a heavy-ion collision at high energy, and to provide a useful tool to 
be applied in numerical calculations of bound state spectra. This was achieved by means of a kinetic 
equation with a collision term consisting of formation and dissociation parts. The equation was solved 
exactly in closed form and the result given in eqs. ()12|) and (|26[). 

The physical content of the solution was discussed and a simplified version of it was given with 
eq. (|31|) in order to illustrate that the yield of bound states form in the medium satisfy the scaling 
embedded in the coalsecence formula, i.e. direct proportionality with the number of bound state 
constituents and inverse proportionality with the system volume. The full solution, on the other 
hand, provides a convenient tool for future numerical studies, which require quantitative knowledge 
of the phase-space distributions of the degrees of freedom involved and of the transition probabilities 
among those degrees of freedom. 
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